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1. INTRODUCTION 

. Tr»r>px will carry a microwave altimeter 

The NASA Ocean Topography Experimen P ? height The capability can be fully 

accurate to a few centimeters for ^detZS to 15 cm or better [Born, et al. 

exploited only if Topex altitude can be indepe Y be tuned with selected nine 

1985], This in turn requires an accurate gravity data type, collected in the first six 

10-day arcs of laser ranging, which will be experimental GPS flight receiver capable 

months of Topex flight. Topex w a s^carry ^ horizon to demonstrate the capability of 

of simultaneously observing six GPS “te determination of the Topex orbit. It has been 

GPS carrier phase and P-code pseudorange P . h mere two -hour arc of GPS tracking 

found that sub-decimeter orbit accuracy can be t ac ^ ^ ^ more ground tra cking sites 

data, provided that simultaneous also valuable for refining the gravity 

[Yunck, et al 1986]. The precision GPS data f rom J° f ^Xitv tuning using GPS measurements, 
model. This paper presents an efficient techmq l ff for f ar fevver gra vity parameters 

Unlike conventional global gravity tuning, - * _ s . $ } it y anomalies which can later be 

in each filter run. These gravity a global gravity map. No 

fo- P such combining. The describes .he approaches 

SSTn 2s study and presents preliminary results of a covanance analysrs. 


2. GRAVITY ANOMALY INFORMATION CONTENT IN TOPEX GPS MEASUREMENTS 

,n this section a theoretica, formulation is derived 

S 2S coordinate (Fig" ->. fi , the radiaf. , 

i r .i. _ j — « ^-T/ni/ ran d€ written 3S 
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Fig. 1 • Rotating h-c-t coordinate 


At, the corresponding uncertainty in the determination of the force 
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Previous analyses have indicated that Topex orbit determination with 
yields 


GPS measurements 


a h / , ^Ah , , <j££ < 0.5 mm/ sec 


Hence, with an altitude h - 1,334 km and nominal velocity C = 7.2 km/sec for Topex, and for a 
samp ing time of At - 2 minutes, local gravity anomaly along a Topex flight path can be determined 

to an accuracy of better than 0.6 mgal in the radial component and better than 0.5 mgal in the other 
two components. 4 * 6 


3. SIMULATION ANALYSIS 


Next, a simulation analysis was performed to numerically estimate the accuracy with which a 
given, but assumed unknown, gravity anomaly can be determined by solving for the coefficients of a 
number of spherical harmonics. In this analysis, a network of six globally distributed tracking sites 
?' eSe inclu< l 0 :nree NASA Dce P Space Tracking (DSN) Sites at Goldstone, California; 

Tu d -! d \S ain ^ n / Ca r berra ' . AuStraha; and three other siles at J a P an ' Brazil ar, d South Africa, 
e ruth model for the gravity anomaly is assumed to consist of a sum of eleven zonal harmonics 

Ji5/ )i6f — / S25r ^ch with a normalized magnitude of 1(H\ j n the filtering process different subsets of 
these terms are estimated and the lumped effects of the selected terms are computed and compared 
with the truth model to assess the accuracy with which the local gravity anomaly can be 
recovered. Other assumptions used in the analysis are included in Table 1. 


Data Type: 

Data Span: 

Data Intervals: 
Gravity Anomaly: 
Station Coordinates: 

Clock Bias: 

Carrier Phase Bias: 
GPS Epoch States: 
Topex Epoch States: 
Zenith Troposphere: 


0.5-cm GPS Carrier Phase 
2 hrs 

2 minutes 

J 15/ Ji6' — / h 5 U0- 6 each, normalized) 
DSN Sites: 3 cm each component (fixed) 
other: 10 cm each component (adjusted) 

3 jisec (adjusted as white noise) 

10 km (adjusted) 

3 m; 0.3 mm/ sec (adjusted) 

10 m; 1 cm/ sec (adjusted) 

2cm 


Table 1. Simulation Analysis Model 

Fig. 2 shows the results when four different subsets of gravity terms are estimated. Note that 

iX Can ** achieved whcn P ro P er model is used in the estimation (Fig. 
2b . °" other hand, a larger discrepancy occurs if the model is not flexible enough to 
. •, e , rUt lg ’ 2c and ^ d ^' ^his leads us to the use of a piecewise constant model for its 
X I tk T prOVlded that the ste P size is fine enough to follow the variation in the truth 
f fK S °i^\ 4 c n ^ US * n ^ 3 P^ ecew ^ se con stant model with a 2-minute step size is shown in Fig. 3. 

With the theoreticaf^rXIion. betWeen ^ a " d ^ ‘ rUth m0dd ( ° 54 mgal) is Cons,stent 


4. THE GRAVITY BIN TECHNIQUE 

Following the above encouraging findings, a more general and efficient technique was 
investigated. This technique solves for "gravity bins", which are 3-D positional deviations on 

lopex, with three parameters at each measurement time point over a period of a few hours. The 
epoch state of Topex and other pertinent parameters are also simultaneously adjusted. These are 
related to Topex current state by the following linearized equation. 
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(a) ALL 1 1 TERMS 


(b) 6 LOWER-ORDER TERMS 








dr dr 

r(,)= Sr 0 r » + S^ v « 


+ 50) 


'f'Z'Jl! 't' hC CUrr<!r " position o( Top «. r » and v„ are posilion and velocity at epoch and 50) is 
the grav, y b,n parametet at time t. Since 5 remains the same lor orbits over VepeaSund^ Vacks 

w 1C will occur every 10 days for Topex, information from repeat orbits can be combined to increase 
the estunahon accuracy. Gravity bins 6 for orbits over differenVground tracks aSTnd^endJTfrom 
one another and can be determined separately. Hence only moderate numto on SeTers^re 

Local sravi,y anoma ' y can 1x1 compu,ed ,rom ,he S0lu,,0ns of 


Ag(t) = 6(t) - ^S(t) 


where g is the nominal gravity field. The global gravity anomaly expressed in terms of spherical 
OanTtmatTo" C °" StrUC,ed ,r ° m * a S— simiL io the^mS 


Fig. 4 shows the result of a 
covariance analysis using a 2.5- 
minute bin size over a 2-hour data 
span which is nearly the Topex orbit 
period. Note that local gravity 
anomaly can be determined to an 
accuracy of 0.2 mgal under most 
circumstances with GPS measure- 
ments over one Topex orbit. This 
accuracy will improve monotonically 
by combining multiple data arcs over 
repeat Topex ground tracks [Ww and 
Yunck, 19881. It is anticipated that 
gravity anomalies of medium wave- 
lengths (1,500 to 3,000 km) can be 
greatly refined with GPS measure- 
ments onboard Topex. 



pravi^v' ^ om P utatlon ‘involved in the transformation of global gravity anomaly from a collection of 
g y ms over the entire globe is well known and can be done without the need for a 

pester. HenCe ' tHe gravity bin techn ique is efficient for global gravity recovery Note that 
grav,ty bins so uhon can be directly applied as a calibration for L effeSs of gravity mism^eltne 

technimie m T detCrmi " at,0n ' CurrentI y under investigation is the application of the gravity bin 
technique to Topex with the baseline (laser ranging) tracking data type, where large data gaps 
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